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Abstract 

We treat the Skyrme model with the breathing mode in a situation 
involving two quartic terms. It is seen that there is a new limit for 
large e due to the breathing mode not found in the usual rotating 
hedgehog. 
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1 Introduction 



Recently many discussions have appeared in the literature [1 -4] about the 
stability of a chiral soliton without the Skyrme term, in particular the soliton 
solution to the nonlinear sigma model by means of the so called breathing 
mode (BM) 0, H| . In this kind of solution, the radial variable, present in the 
hedgehog solution, is scaled by a dynamical variable so that its subsequent 
quantization furnishes an extra energy for the system. The proposal is that 
this energy is responsible for stability in a mechanism similar to the hydrogen 
atom. In fact, because of the BM, one can see in a rough approximation, 
using the uncertainty principle, that the classical energy (and radius), even 
when the rotations are absent, has a minimum. Although the results are 
good and comparable to that of the Skyrme model f| [| (which are obtained 
when as well as the nonlinear sigma term, the Skyrme term is also present) 
there are some problems that remain, one of them being the choice of the 
profile function F(r) (determined in the Skyrme model by the minimization 
of the static mass). When the Skyrme term is absent, the static mass has no 
minimum and the profile function is quasi-arbitrarily chosen. This fact can 
lead to a new kind of instability as one can choose a function that results in 
divergent integrals even when the required boundary conditions are obeyed 
[0], so that the presence or not of the BM stabilizer term would make no 
difference. 

As the BM has shown its importance in the nonlinear sigma model through 
the existence of a special stable soliton configuration, one can argue what is 
its relevance in the Skyrme model (for example, how accurately one can de- 
scribe the nucleon and delta resonances). Although there are already some 
treatments of this subject in the literature || [|, mainly about the stabil- 
ity [|Kj] of the Skyrme model with a vibrational solution, here we intend to 
explore the relationship between the Skyrme model and the nonlinear sigma 
model, more specifically how one can get the results of the nonlinear sigma 
model with BM as was done, for example, by Carlson M, from the resolu- 
tion of the Schrodinger equation of the Skyrme model with BM. In other 
words, we depart from the Skyrme model, which has a well defined profile 
function, perform the quantization of the rotational and vibrational degrees 
of freedom and take a suitable limit on the subsequent Schrodinger equation. 
As a result, we get the vibrating nonlinear sigma model with finite integrals. 
This kind of limit is a property only of the vibrating skyrmions. The same 
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limit exists for all static properties. This is the main idea of the present 
contribution. 

2 The breathing mode 

The SU(2) Skyrme model in a version that includes the two independent 
fourth order terms (Skyrme and symmetric quartic terms (SQT)) and in 
the chiral symmetry limit is described by the following Lagrangian density 

(h = c = iy. 

C sk = -^Tr(L^) + ^ 2 Tr([L„ L v f) + ^{[Tr(L^) 2 ]}, (1) 

with L M = Wd^U, F n =186 MeV the pion decay constant, e and 7 parameters 
to be adjusted and U = e lT ' n( - x ' tS) the pion field. Skyrme introduced soliton 
solutions in his model through the field configuration called hedgehog. In this 
ansatz, the pion field has a radial direction and depends on a function F(r) 
called the chiral angle. This function is determined by the minimization of 
the static mass associated with the hedgehog solution. Well defined states of 
spin and isospin are obtained through a rotation of the classical solution. The 
BM can be taken into account by means of a scaling on the radial variable r: 

U = A(t)e W4F{r > R{t)) A\t) } (2) 

where R(t) is the dynamical variable that characterizes the soliton breathing 
mode and A(t) G SU(2) is the rotation matrix. The use of the solution (Q) 
in the Lagrangian density (|1|) and an integration over all space leads to the 
following Lagrangian: 

L = R 3 bd^ + R(k - 7 m)dX - Rc - - ~ 7 * + R 2 (f - 771) + RR 2 d, (3) 

R 

with a M the four variables that parametrize the rotation matrix A(t) and the 
constants are listed in appendix. We stress that the integrals c, b and d are 
the contributions from the nonlinear sigma term to the static mass, to the 
rotational term and to the BM, respectively. In the same order, h, k and / 
are the contributions from the Skyrme term and Z, m and n the contributions 
from the symmetric quartic term. 
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One can pass to the Hamiltonian formulation through a Legendre trans- 
formation, whose result is: 



H 



7T 



+ 



PR 



AR(R 2 b + k — 7m) A(Rd + f - jn) 



+ Rc + 



h-jl 
R ' 



(4) 



with 7r 2 = TT/j,^ the square of the momentum conjugate to a M and pr the 
momentum conjugate to R. Now, we will write the Hamiltonian (f|) through 
a set of reduced variables obtained when we make use of the constraint on the 
rotational variables (a^a M = 1). We remind the reader that if we first make 
the reduction of variables at the classical level and perform the quantization 
after, or first perform the quantization and afterwards restrict adequately 
the vectors of the Hilbert space, we obtain the same result [11], Q . Another 



way is to construct a Hamiltonian compatible with the constraint by means 
of Dirac's procedure for constrained systems In this case there appears 
a contribution to the zero point energy inversely proportional to the moment 
This contribution will not be treated in the present work. 

a\R) and p = {tt^Pr), i— 1,2,3, we 



of inertia [12, 14 



Using the reduced coordinates q 
write the Hamiltonian as: 



H = lp a g ab p b + cq±-' H 1 



T 



where: 



-,ab 



R(R 2 b + k — 7m) 








1 

2{Rd + / - 771) 



(5) 



(6) 



and 



(7) 



We write the quantum Hamiltonian via the DeWitt [15| prescription: 
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or 



J_ f_#J_ 4B 4 + (5(/ - 7w)/d + (fc - ~fm)/b)R 2 + 2(f - jn)(k - jm)/bd d \ 
~4d\dR 2 + R{R 2 + {f- -fn)/d){R 2 + {k- 7m) /6) aRj ^ 



i2 \R(bR 2 + k - 7m) # / # 

where g is the determinant of the metric tensor. We make a transformation 
on the wave function ip as indicated in 0, 

i, = { 2(dfi2+ /~ 7,i) } 1/4 {2fl(^ + k - 7 m)}- 3/ V (10) 

and eliminate the first derivative from the equation. So: 

1 <9 2 
AddR 2 ^ 

8R e + (18(f --fn)/d+10(k--fm)/b)R 4 
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+ Wd ((/ - 771) /d + # 2 ) 2 {{k - 7 m)/6 + # 2 ) 2 ' 

(15(/ - 7n) 2 /^ 2 + 18(/ - 77i)(fc - ^m)/db)R 
+ 16d ((/ -^n)/d+R 2 ) 2 ((k-^m)/b + R 2 ) 2 * 
(k - -fm) 2 R 2 /b 2 + 18{k - 7m) (/ - jn) 2 /bd 2 - 6{f - jn)(k - jm) 2 /db 2 
16d ((/ - 771) /d + # 2 ) 2 ((fc - 7m) /& + i? 2 ) 2 ' 
(i? + (/- 7 n)/d) / + fr- 7 A 

+ i2 \R(bR 2 + k - 7m) # ^ 

_ i? 2 + (/- 7 n)/d 



(11) 



In order to obtain a solution to this equation, we observe that the behavior 
of near the origin is: 



cj)(R^0)~R* I* a "-7™ yj (12) 

and we require that 4> goes to zero at infinity. 
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3 The limit of large e 



Usually in the Skyrme model, the two free parameters of the theory are 
chosen to reproduce the nucleon and delta masses. As a result, the value 
of e that reproduces such masses does not correspond to the e value that 
minimizes them. This is not so important in the context of the value of e to 
be used, because there is no necessity to use the e that minimizes the mass. 
But it is important in the sense that it shows that for e going to zero or e 
going to infinity, the mass diverges (see figure [I]). Now, the question is what 
happens with the mass in that limit when along with the rotations, the BM 
is also present? The answer is that for small e the mass still diverges, but 
for sufficiently large e, the mass becomes a constant in e and all works as if 
only the nonlinear sigma term were present, with a profile function F that 
minimizes the static mass of the Skyrmion. 

The results for the nucleon and delta masses and their excitations are the 
same as those obtained some time ago by Carlson . In that work, Carlson 
pointed out that it is possible to obtain stable soliton configurations from the 
nonlinear sigma model when the Skyrme term is used as a constraint in the 
action of the system. This was shown using scale invariant arguments in the 
functional integral for the field configurations U, dividing the configurations 
into equivalence classes and extending the functional integral to the equiva- 
lence classes of configurations. The resulting functional integral depends of 
a dilatated field and involves a restriction on a particular function G(U) of 
U. This function must preserve the chiral limit and must not be invariant 
by dilatations as, in this case, it would not be possible to choose a field con- 
figuration that was representative of its class. Carlson chose for G(U) the 
Skyrme term. With this choice, it is possible to solve the Euler-Lagrange 
equation without the explicit presence of the Skyrme term in the action. So, 
Carlson obtains a hamiltonian of the nonlinear sigma model with integrals 
dependent of the chiral angle from the Skyrme term. 

Now, we will call y the dimensionless variable that scales the radial vari- 
able in the nonlinear sigma model. It is related to R by R = ey. If one makes 
this replacement in the Schrodinger equation @ is possible to verify that for 
6 — ► oo, it reduces to: 

4yd { dy 2 y dy } ^ by 3 ^ c ^ (-^) 
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which is e independent. As the SQT and the Skyrme term have the same 
e dependence, the equation has the same functional form when 7 is zero or 
not in the limit of large e, with the only difference being the function F 
that minimizes the classical mass. The expression (|T3| ) is the Schrodinger 
equation for the BM when only the nonlinear sigma term is present. In fact, 
the numerical resolution of equation (|9]) shown in figure ([I]), tells us that for 
e > 80 the mass is already that calculated by Carlson. The other static 
properties have the same kind of limit. As an example, we write down the 
expression for the axial coupling g a . It is easy to verify that g a , given in the 
vibrating Skyrme model by 



1 w fl? V VSdR f 2 sm(2F) \ 

87T f°° , . 

+ — / dpsin 2 FF' 
6e 2 Jo 

+ ( r dppsin{2F){F'f + r dpsin 2 F S ^^\ , (14) 
Se 2 \Jo Jo p J 

is reduced, in the limit of large e, to: 

7T / y 2 I ip | 2 ^/gdy 2 / sin(2F)\ 



9a 



3 7Wv5d» y. dpp r (15) 



which is the corresponding expression in the nonlinear sigma model. The 
same holds for all other static properties. 



4 Numerical results 

Figure (HD shows the mass curve for the case of the rotating hedgehog with 
the Skyrme term (case I), Skyrme term plus BM (case II) and Skyrme and 
SQT terms plus BM (case III). The value used for the pion coupling constant 
is Fn = 186 MeV and 7 is 0.11 which is the maximum 7 value that mimimizes 



the classical mass. This value agrees with those found in the literature |px 
The discrepancy between cases II and III is due to the presence or absence 
of the SQT. As seen in equation (|TTD, the SQT contributes negatively to the 



energy so that for low e, the mass calculated when the SQT is present is 
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lower than that calculated when only the Skyrme term is present. In the 
high e region, the only difference is in the constants b, c and d of equation 
( P~5| ) which are calculated with different profile functions. For case II they 
are 125.01/F w e 3 , 18.23F^/e and 160.79/i^e 3 and for case III 52.57/i^e 3 , 
U.72F w /e and 72.87/i^e 3 respectively (see table in the appendix). The 
asymptotic mass value for case II is 1101 MeV and for case III 1152 MeV 
compared with an experimental value given by 939 MeV. 

As the energy is not absolutely defined in our treatment, we decided to 
search for mass differences (fig. 0). Initially we have two free parameters in 
the model. One of them, as stated before, was chosen to fit the pion decay 
constant. The other, in principle, could be chosen to fit the mass diference 
between the states with j = 1 / 2 and j=3/2. The problem is that for small e 
the mass difference is too small and for large e it is independent of e so that 
in this region we have only one free parameter. As we choose to fit F n , the 
biggest mass difference for case II is 285 MeV and for case III is 306 MeV. 
The experimental value is 293 MeV. The differences between cases II and III 
arise because the centrifugal term in III is bigger than II implying that the 
mass difference is also bigger. In figure (0) we also plot the mass difference 
for the rotating hedgehog with the Skyrme term. As is well known, it behaves 
like e 3 . Finally, to illustrate the behaviour of the other static properties, the 
isoscalar radius curve is plotted in figure (^). One can see that for large e, the 
isoscalar radius for the rotating hedghog goes to zero, as expected, because 
in that region the effect of the Skyrme term is weak and the soliton collapses. 
On the other side, when the BM is present, it is verified that for large e the 
isoscalar radius goes to a constant value. For case II, the mean charge radius 
is 0.56 fm and for case III is 0.58 fm. Another example is the isoscalar mean 
magnetic radius. For case II, it is 0.96 fm and for case III is 0.92 fm. 

5 Conclusions 

We have studied the Skyrme model in a solution for the pion field that in- 
cludes both rotational and vibrational degrees of freedom. It was found that 
the presence of the breathing mode in the solution allows one to take the limit 
e — > oo - a situation not found in the usual rotating hedgehog. As a conse- 
quence, we recovered the earlier results of Carlson p|. We did this study in 
two situations: one involving only the Skyrme term and the other involving 
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the Skyrme term plus the symmetric quartic term. The limit for large e is 
the same for the two cases with the only difference being the function F that 
minimizes the classical mass. Some static properties were presented and it 
was noticed that the masses of the nucleon and delta states as well as the 
mass difference between them are larger when the symmetric quartic term is 
present. This is due to centrifugal effects. We also think that the method 
presented here provides some justification for the use of a particular chiral 
angle in the model proposed by P. Jain et al. jl). 



One of us (F.M.S.) acknowledges helpful discussions with M. Betz and 
G. L. Thomas. This work was supported by the Brazilian agencies CNPq, 
FAPERGS and FINEP. 

Appendix: Integrals containing the chiral 
angle 

Let F(p) be the chiral angle, s = sinF, c = cosF, F^ the pion decay constant 
(71 — > p + v) and e the Skyrme dimensionless parameter. So, the explicit 
form for the constants appearing in the Lagrangean is given by: 




(1) 









9 = 



Se 3 ^ Jo 



16n r 



00 
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dF 
dp 



) 



2 



+ -» P 2 dp, 




(5) 
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Integrals 


7 = 


7 = 0.11 


b 


125.0145 


52.573 


c 


18.23082 


14.716 


d 


160.7893 


72.874 


f 


43.91 


39.994 


g 


88.52515 


75.37 


h 


18.23 


22.965 


1 




74.988 


m 




194.834 


n 




193.077 



Table 1: Numerical values for the integrals containing the chiral angle when 
only the Skyrme term is present (7=0) and when the SQT (7=0.11) is added. 
The parameters e and F n are not included. 




The coefficients c, h and / are the contributions for the soliton static mass 
from the nonlinear sigma term, the Skyrme term and SQT respectively. In the 
same way, the coefficients b, g and m are the contributions for the rotational 
energy and d, f and n the contributions for the vibrational energy. 
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Figure 1: Energy curve for the rotating hedgehog with and without BM 
in function of the Skyrme parameter e. In cases II (dashed line) and III 
(continuous line), where the BM is present, the energy decreases with e until 
a region where it becomes e independent. In case I (dotted line), the energy 
diverges as e tends to zero or infinity. 
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Figure 2: Mass difference curve between the states j—1/2 and j—3/2 for the 
rotating hedgehog with and without BM in function of the Skyrme parameter 
e. In cases II (dashed line) and III (continuous line), the mass difference 
increases until a region where it becomes e independent. In case I (dotted 
line), the mass difference diverges as e goes to the infinity. 
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Figure 3: Isoscalar radius curve for the rotating hedgehog with and without 
breathing mode in function of the Skyrme parameter e. The notation is the 
same as in figure 1. In the case that the breathing mode is present, the 
isoscalar radius decreases with e until it reaches a minimum and tends to 
a constant. In the other case, the radius goes continuously to zero as the 
strength of the Skyrme term becomes weaker. 
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